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Small Numerators Canceling Small Denominators: Is Dyson's Hierarchical Model 

Solvable? 

Y. Meurice 

Department of Physics and Astronomy, The University of Iowa, Iowa City, Iowa 5224-2, USA 

We present an analytical method to solve Dyson's hierarchical model, involving the scaling vari- 
ables near the high-temperature fixed point. The procedure seems plagued by the presence of 
small denominators as in perturbative expansions near integrable systems in Hamiltonian mechan- 
ics. However, in all cases considered, a zero denominator always comes with a zero numerator. 
We conjecture that these cancellations occur in general, suggesting that the model has remarkable 
features reminiscent of the integrable systems. 
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In many physical problems involving nonlinear flows, 
a common strategy consists in constructing a system of 
coordinates where the flow becomes linear. The action- 
angle variables in Hamiltonian mechanics provide well- 
know examples of such a procedure. Whenever the angle 
variables can be constructed, they evolve linearly with 
time and expressing the original variables in terms of the 
new ones solves completely the original problem. The 
problems for which well-defined angle-action variables 
can be constructed (e.g. Kepler's problem or the free 
rigid body) are very distinguished and called integrable 
systems. A large number of numerical experiment has 
lead us to believe that in a generic way, small perturba- 
tions destroy integrability. This point of view was was 
first inferred by H. Poincare who pointed out the exis- 
tence of small denominators in the canonical transfor- 
mation designed to eliminate the angle dependence of a 
perturbed Hamiltonian. 

In this letter, we discuss the question of small denom- 
inators for renormalization group (RG) flows. The vari- 
ables which play the role of angle variables are the scal- 
ing variables introduced by Wegner Q. Near a fixed 
point, the RG flows can be linearized. The problem of 
expressing the physical quantities in terms of variables 
which transform as in the linear approximation when the 
non-linear terms are taken into account, is analogous to 
removing the angle dependence of a perturbed Hamilto- 
nian. If the task can be carried through, one obtains 
analytical expressions for the RG flows. However, as we 
will show, small denominators appear. Does this mean 
that, as in Hamiltonian mechanics, in generic situations 
the construction will fail? 

Surprisingly, we found in a numerical calculation per- 
formed with Dyson's hierarchical model [^|J^], that zero 
denominators were systematically canceled by zero nu- 
merators. These remarkable cancellations suggest that 
either the model considered is as distinguished as the in- 
tegrable systems of classical mechanics or that there ex- 
ists a general mechanism that allows us to circumvent the 
small denominator problem for RG flows. The example 
of the two-dimensional Ising model shows the importance 



of having a non-trivial model which can solved in closed 
form. The results presented below indicate that Dyson's 
hierarchical can be solved analytically. 

During the last decades, the RG method has been suc- 
cessfully applied to many important problems in field 
theory and statistical mechanics: the critical behavior 
of ferromagnets and superconductors, the confinement 
of quarks or the generation of mass for the W and Z 
bosons. However, its practical implementation is still 
a formidable enterprise. Typically, expansions are of- 
ten available near fixed points, but not to orders large 
enough to allow one to extrapolate between fixed points. 
Unfortunately, the calculation of physical quantities (e. 
g. the magnetic susceptibility) beyond an order of mag- 
nitude estimation, requires a calculation of the flows in 
crossover regions. One has then to rely to Monte Carlo 
simulations to achieve this goal. One also needs to select 
a small set of interactions which closes reasonably well 
under RG transformation near both fixed points. Inter- 
esting examples of such lattice Monte Carlo calculations 
are given in Refs. Q for scalar theory and || for gauge 
theories. 

In order to get analytical results, further approxima- 
tions are needed. One possibility consists in using hierar- 
chical approximations such as the one derived by Wilson 
in Ref. || and resulted into the "approximate recursion 
formula" . In this approximation, only the local inter- 
actions get renormalized and the flow can be calculated 
from a simple integral formula. Retracing Wilson's con- 
struction from the beginning, one can restore the other 
renormalizations perturbatively. In order to perform this 
task, one would like to have a closed form solution in the 
hierarchical approximation. In order to achieve this goal, 
we have proposed to use the Fourier representation of 
the integral formula and to construct the scaling vari- 
ables in this basis. In this process, we identified the ex- 
istence of small denominators possibly ruining the whole 
approach. This problem can be avoided in special cir- 
cumstances, for instance for flows starting exactly along 
the unstable direction of a non-trivial fixed point. With 
this restriction, we found M expansions with overlapping 



1 



domains of convergence in the crossover region. 

In the following, we discuss the problem of small de- 
nominators in the construction of the scaling variables 
near the high-temperature (HT) fixed point of Dyson's 
hierarchical model. The treatment of the this model is 
mathematically similar to the one of the approximate re- 
cursion formula. However there exists a large literature 
on Dyson's model [0-1C], the non-trivial fixed point is 
known very precisely 1 1 1 1 and the HT expansion studied 



to a very large order [ 12] . 

For a description of this model as a spin model, we 
refer to Ref. 13 13], while the details of the derivation 
of the RG flows as expressed below can be found in . 
To make a long story short, all the information regard- 
ing the local interactions after n RG transformations is 
encoded in a function R n (k) = 1 + a n .ik 2 + a n ,2k A + . . .. 
The logarithm of this function generates the connected 
zero-momentum Green's functions at finite volume. The 
recursion formula reads 



R n +i(k) = C„+i exp 
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We fix the normalization constant C n so that i?„(0) = 1. 
We use the parametrization c = 2 1 ~ 2 / D which implies 
that a free massless field scales in the same way as in 
a usual D-dimensional theory. In this formulation, the 
temperature dependence has been absorbed in the initial 
Ro(k). For an Ising measure, Ro(k) = cos(\f]3k), while 
in general, wc have to numerically integrate the Fourier 
transform of the local measure to determine the coeffi- 
cients of Ro(k) expanded in terms of k. In general, a ni i 
is of order f5 l in the HT expansion. 

In the HT phase, polynomial truncations of order l max 
in k 2 provide rapidly converging approximations [p"3| . 
The RG flows can be expressed in terms of the quadratic 
map 
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U n ,l 
U n ,0 



with 



and 
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{n + v-a)\{2a)\ 



(2) 



(3) 



(4) 



for fi + v > a and zero otherwise. We use "relativistic" 
notations. Repeated indices mean summation. The greek 
indices [i and v go from to l ma x, while latin indices i, 
j go from 1 to 

The diagonalization of the linear RG transformation 
near the HT fixed point is quite simple because it is of 
the upper triangular form. From Eq. (|J), one finds the 
spectrum 



A (r) = 2(c/AY 



(5) 



in agreement with Ref. 



Using the matrix of right 



eigenvectors, M.\i>\ = Xmi/jJ, we introduce new coor- 
dinates such that a nt i — tpfh^r. This diagonalizes the 
linear RG transformation. Note that the form of the 
eigenvectors guarantees that h n> i is also of order fr . More 
details are given Ref. , where all the quantities related 
to the HT fixed point are dressed with a "tilde" omitted 
here . In summary, the RG flows in the new coordinates 
can be written as 



Ki+i.i — 



1 + APh ri 



^0 h"n,phn,q 



(6) 



with coefficients calculable from Eq. (Q). 

We now express the h n j in terms of the scaling vari- 
ables y n ,i, . . . , y n ,i maa: which transform under a RG trans- 
formation as y n +i,i — KUn.i- If we can construct func- 
tions hi and yi such that h n j = hi(y n ) and y n ,i = Jft(h-n), 
then we get a complete analytical expression of h n j 
(which contains all the thermodynamical quantities) in 
terms of hoj (which depends on the initial energy den- 
sity): 



K,X = M A "yi(ho), A 2J/2(ho), . . .) 



(7) 



The feasibility of this approach is demonstrated for a 
one-dimensional example in Ref. Jl^] . 

We now discuss the construction of the hi . We use the 
expansion 



hi = yi 



Sl,ixi2-VlV2 



(8) 



where the sums over the i's run from to infinity in each 
variable with at least two non-zero indices. In the follow- 
ing, we use the notation i for . . .). Plugging the 
expansion into Eq. and requiring that the one step 
advance is obtained by rescaling the scaling variables by 
their associated eigenvalue, we obtain 



si 



(Tim A (m) ~ A ' 



(9) 



with 



= Ej+ k =i(Af%J^k ' S U Urn Afc)A%, k ) (10) 
— Sj+k+r=i S; J Tim ^(m) ^0 ?S P,k s 9,r ■ 

For a given set of indices i, we introduce the notation 

l q {i) = Y J i m m q . (11) 

TO 

One sees that Iq is the degree of the associated monomial 
and X\ its order in the HT expansion (since yi is also of 
order /3 ) . Given that all the indices are positive and that 
at least one index is not zero, one can see that if j + k = i 
then 2q(j) < T q {\) and I q (k) < 2"<j(i)- Consequently, Eq. 
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( |lO| ) provides a solution order by order in Xq or in T\ 
(since the r.h.s is always of lower order) assuming that 
the small denominator problem can be avoided. 

Using the parametrization c = 2 1 ~ 2 /- D , the zero de- 
nominators in Eq. (H) appear when 



D~l{D + 2) = DJ Q -(D + 2)X x 



(12) 



If D and D + 2 have no common factors, this equa- 
tion has non-trivial solutions for indices such that To = 
(2 + D)q + 1, for q a strictly positive integer and I such 
that Zi = Dq + I. If D and D + 2 have common factors, 
we can proceed in the same way but with these common 
factors removed from both D and D + 2. 

We have investigated numerically the well-studied case 
D = 3 with a polynomial truncation l max = 25. For prac- 
tical reasons, we have limited our study the case where 
hi depends only on the two leading variables j/i and j/2- 
This restriction is self-consistent since if we start for in- 
stance with 2/3 = 0, the multiplicative renormalization 
of the scaling variables implies that this condition stays 
valid after n iterations. Since — = . . . = 0, we 
will use notations such as s;^^ or /^(j/i, 2/2)- In ad- 
dition, we have concentrated our attention to the large 
order behavior in the leading variable 2/1 and calculated 
the first order corrections in the subleading variable 2/2 If 
i2 = 0, we have the following sets of small denominators: 
i\ = 5q + 1 and I = 2q + 1. If i% = 1, we have i\ = bq 
and I — 2q + 2. We have calculated the numerators corre- 
sponding to these zero denominators for q = 1, 2, 3 and 
4. 

We have checked our calculation of the hi with two 
different methods. First we have considered a configura- 
tion h n ^i = hi (2/1,2/2) for particular values of 2/1 and 2/2 
and then calculated the one step backward configuration 
h n ~\,i = /iz(A^f 1 yi, A 2 ~ 1 2/2)- We then calculated the one 
step forward using these h n -\ t i and the exact formula 
Eq. (^|). Comparison between the two h n ^i for values of 
2/1 and 2/2 varying between 0.1 and 0.001 and for various 
truncations in the power of 2/1 and 2/2 considered, showed 
errors scaling like the powers neglected with coefficients 
compatible with the order of magnitude of the coefficients 
involved in the expansion. Second, we used the expres- 
sion yi(hi) calculated in Ref. up to order 11 in [3 and 
checked J/i (/iz (2/1 , 2/2 ) ) = Vi with errors ranging between 
10 -14 and 10~" 16 for the various coefficients of the higher 
order terms. 

In Fig. 1 we show the absolute value of N% t jfi. The 
calculations have been performed with different arith- 
metic precisions. We have used Mathematica 4.0 with 
$MaxPrecision and $MinPrecision both set to the 
same value pr. We have considered the cases pr=20, 
30 and 40. The graphs shows A^.o is more than twenty 
orders of magnitudes smaller than its peers when pr = 20 
and then drops by ten orders of magnitude each time the 
precision is increased by 10, while all the other values 
stay stable. This is a strong evidence for ^3^,0 = 
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FIG. 1. Logio(\N 3 jfi\) versus j. 

In Figs. |^ and|^, we show Ni il and Nij u i calculated 
with pr=20. One sees that A^+i, 5<?+i,o and -^2,3+2,5^,1 
are more than 20 orders of magnitude smaller than the 
naive interpolation. We have checked in each of these 
cases that the small value drops by ten order of magni- 
tudes each time the precision is increased just as in Fig. 

0- 



O) 

o 




10 15 20 25 

j 

FIG. 2. LogioQNij l0 \) versus j for / = 1, . . . 9. 
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FIG. 3. Log 10 (\Ni JA \) versus j for I = 1, ... 10. 



3 



In summary, for each zero denominator considered, we 
found a zero numerator. We thus conjecture that all nu- 
merators corresponding to zero denominators are zero. If 
this conjecture is correct, the corresponding to non- 
zero denominators are calculable from Eq. (|l0|) and those 
corresponding to zero denominators are undetermined. 
The calculations performed above have been done with 
these coefficients treated as undetermined. The figures 
have been drawn with these coefficients set to zero. This 
choice is not essential, we have considered other choices 
where the undetermined coefficients have been set to val- 
ues of the same order as the coefficients above and below 
(in «i) and reached identical conclusions. This is a non- 
trivial statement. For instance, the undetermined coef- 
ficient S3. 6,0 appears explicitly in the numerator of the 
equation for 35,11,0, however it is multiplied by a very 
small number which drops when pr is increased. 

We have checked that the individual terms in the zero 
numerators were not zero. The fact that we were able to 
obtain very precise cancellations without any fine-tuning 
suggests the existence of closed form formulas or of a 
symmetry forbidding these terms. More generally, the 
equally spaced spectrum (on a log scale), the number of 
terms at a given level m (the order in HT expansion) 
equal to the number of partitions of m Q and the exis- 
tence of "nested" ambiguities are somehow reminiscent 
of string theory. A possible starting point to discover 
these hypothetical symmetries would be to exploit the 
fact that since for instance 7/3 and yf transform the same 
way under a RG transformation, there exist ambiguities 
in the construction of hi in terms of the scaling variables. 

We should say a few words about the small denom- 
inators near other fixed points. The spectrum of the 
linearized RG transformation near the non-trivial fixed 
point [|ll| can be calculated numerically |l(J . Investigat- 
ing the small denominators of the form Xf ~ \ r m for the 
first twenty eigenvalues and for powers not larger than 
100. The best solution we found was A2 — A4 with two 
parts in a thousand. A more detailed study is necessary 
to decide if the rate of decay of the coefficients is sufficient 
to take care of the smaller denominator which will appear 
at larger orders. On the other hand, the spectrum at the 
Gaussian fixed point Q is Xj = 2c~ : > . There are many 
zero denominators in integer dimensions, e.g., Ai = A| 
for D — 3. If the numerators are not zero, one can "re- 
pair" jl] the situation by considering n-dependent coef- 
ficients. This excludes a solution of the form of Eq. (Q) 
but it generates logarithmic corrections which are neces- 
sary. These can even be observed in the large order of 



the HT temperature expansion in Ref. |17| . 

In conclusion, we have found remarkable cancellations 
of small denominators by small numerators. If these can- 
cellations occur in general, it is possible to calculate an- 
alytically the thermodynamical quantities of Dyson's hi- 
erarchical model in the HT phase using Eq. (|l0|). Our 
results suggest that this model has features analogous 
to the integrable systems in Hamiltonian mechanics. As 
such, the model would stand out as a first approximation 
to be used in situations where the conventional pertur- 
bative expansions are not reliable. 
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